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ABSTRACT 

We present the results of the northern extension of the HI Parkes Zone of Avoidance Survey, a blind HI 
survey utilizing the multibeam receiver on the Parkes 64-m telescope. In the two regions studied here, I = 36° 
to 52° and I = 196° to 212°, \b\ < 5°, we have detected 77 HI galaxies, twenty of which have been previously 
detected in HI. The survey has a median rms noise of 6.0 mJy beam" 1 and is complete to a mean flux density 
of 22 mJy. We have searched for multiwavelength counterparts to the 77 galaxies detected here: 19, 27, and 1 1 
have a likely optical, 2MASS, and IRAS cataloged counterpart, respectively. A further 16 galaxies have likely 
visible counterparts on the Digitized Sky Survey. The detection of these 77 galaxies allows a closer inspection 
of the large-scale structures in these regions. We see several filaments crossing the Galactic plane, one of which 
appears to be the continuation of a sine-wave like feature that can be traced across the whole southern sky. An 
analysis of the HI mass function suggests that the regions studied here may be underdense. One particularly 
noteworthy galaxy is HIZOA J0630+08 (l,b = 203°, -0.9°) with a velocity of 367 ± 1 km s _1 . We suggest that 
it belongs to the nearby Orion Group which includes a small number of dwarf galaxies. The newly detected 
galaxies improve our understanding of the properties of several voids, such as the Orion, Gemini, and Canis 
Major Voids. 

Subject headings: galaxies: distances and redshifts — galaxies: fundamental parameters — large-scale struc- 
ture of universe — surveys 



1. INTRODUCTION 

The obscuring effects of dust and high stellar density in the 
Milky Way have historically prevented the detection of galax- 
ies that lie behind the disk of our Galaxy. This zone of avoid- 
ance (ZOA) typically comprises more than 20% of the optical 
extragalactic sky and nearly 10% of the infrared extragalactic 
sky, although recent surveys have begun to lessen its extent. 
The presence of the ZOA does not affect extragalactic popu- 
lation studies, as galaxies in the highly obscured ZOA should 
not differ from those detected in regions of low obscuration. 
A complete census of galaxies across the Galactic plane is 
necessary, however, to fully understand the dynamics of the 
Local Group as well as the large-scale structure of the Local 
Universe. 

Many early ZOA studies focused on finding a highly ob- 
scured, nearby, massive galaxy. Several such galaxies have 
been found, including Maffei 1 and 2 (Maffei 1968), Circinus 
(Freeman et al. 1977), and Dwingeloo 1 (Kraan-Korteweg 
et al. 1994). While the IC 342/Maffei group and other 
nearby massive galaxies largely define the motion of the Lo- 
cal Group, the direction of the Local Group acceleration is 
strongly affected by the full ZOA mass distribution out to 
6000 km s" 1 (Kolatt, Dekel, & Lahav 1995), and possibly be- 
yond. 
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Several attempts have been made in recent years to sam- 
ple the ZOA more completely. Optical galaxy candidates 
have been identified by searches of the Palomar Observatory 
Sky Survey (POSS) prints and Schmidt atlases (see Kraan- 
Korteweg & Lahav 2000 for a review). Highly obscured 
galaxies have also been detected in the near-infrared (NIR) 
and far-infrared by the Deep Near-Infrared Survey of the 
Southern Sky (DENIS; Schroder, Kraan-Korteweg, & Mamon 
1999), the Two Micron All Sky Survey (2MASS; Jarrett et 
al. 2000), and the Infrared Astronomical Satellite Sky Survey 
(IRAS, e.g. Nakanishi et al. 1997). 

While useful in reducing the size of the ZOA, optical and 
infrared surveys are ultimately limited by Galactic extinction 
and source confusion. The most thorough optical searches are 
complete to an apparent diameter of D = 14" and a magnitude 
of B 2 5 = 18 m .5 at extinctions of A B < 3 m .O (Kraan-Korteweg 
2000, Woudt & Kraan-Korteweg 2001). The DIRBE///M5 
100 fim extinction maps (Schlegel, Finkbeiner, & Davis 1998) 
directly measure the dust column density and indicate that 
A B > 3 m .O for 9.5% of the sky, although the data in the Galac- 
tic plane are not yet well-calibrated. The completeness of in- 
frared surveys of the ZOA is affected primarily by source con- 
fusion. Results from the DENIS survey indicate that NIR sur- 
veys can easily detect galaxies down to Ag ~ 10 mag, and the 
brighter ones far beyond (Schroder et al. 1999). The 2MASS 
survey is expected to be complete to 13 m .5 at 2.2 fj,m for the 
majority of the sky, but for 5° < \b\ < 20°, this completeness 
drops to 12 m .l (Jarrett et al. 2000). 

HI emission is not affected by the dust in the Milky Way. 
Consequently, blind HI surveys provide a unique means by 
which to identify gas-rich ZOA galaxies. The first blind ZOA 
HI survey was conducted by Henning (1992) using the Green 
Bank 300-ft telescope. More recent blind HI surveys include 
the Arecibo Dual Beam Survey (ADBS, Rosenberg & Schnei- 
der 2000) and the Arecibo HI Strip Survey (AHISS, Zwaan et 
al. 1997), which intersect with but do not focus on the ZOA. 
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The first systematic HI survey to focus on the optically most 
opaque part of the ZOA is the Dwingeloo Obscured Galaxies 
Survey (Henning et al. 1998, Rivers 2000) which sampled the 
northern ZOA (30° < I < 220°, \b\ < 5°.25) to velocities of 
4000 km s" 1 in the Local Group standard of rest with an rms 
of 40 mJy beam" 1 . Each of these surveys has contributed to 
our understanding of gas-rich galaxies and the structure of the 
local Universe, but all have been limited in either sensitivity 
or sky coverage. 

The 2 1 -cm multibeam receiver (Staveley-Smith et al. 1 996) 
on the Parkes 6 64-m telescope is the first instrument capable 
of conducting both sensitive and large-area blind HI surveys 
in a reasonable time. Two multibeam surveys commenced in 
1997. The HI Parkes All-Sky Survey (HIPASS, see Meyer et 
al. 2004) is a blind HI survey of the sky south of 6 = 25° and 
has an rms noise of 13 mJy beam" 1 . The HI Parkes Zone of 
Avoidance Survey is a more sensitive (rms = 6 mJy beam" 1 ) 
survey of the region of the Milky Way accessible from Parkes 
(/ = 196° to 52°, \b\ < 5°). An analysis of shallower data for 
the southern ZOA (rms = 15 mJy beam" 1 , 1 = 212° to 36°) led 
to the discovery of 1 10 galaxies, 67 of which were previously 
unknown (Henning et al. 2000). Here we report on the results 
for the complementary region north of Dec. 0° (/ = 36° to 52° 
and / = 196° to 212°). The results reported here are, however, 
at full survey sensitivity (rms = 6 mJy beam" 1 ). Although this 
region overlaps with the Dwingeloo survey, the Parkes survey 
is much more sensitive and represents the first probe, beyond 
a few Mpc, of this region of the Universe, which lies near both 
the Local and Microscopium Voids. Data of similar depth for 
the remainder of the ZOA (/ = 212° to 36°) will be reported 
elsewhere. 

In §2, we describe the observations and the method by 
which the data were reduced. The galaxy search method, 
measurement of HI parameters, search results, and survey 
completeness are discussed in §3. In §4, we describe the re- 
sults of a search for multiwavelength counterparts to the sur- 
vey galaxies. The large-scale structures revealed by the sur- 
vey are discussed in §5, and the mass function of the galax- 
ies is discussed in §6. Throughout the paper, we assume 
Ho = 75 km s" 1 Mpc" 1 . 

2. OBSERVATIONS AND DATA REDUCTION 

The observations for the northern extension of the HI 
Parkes ZOA survey began in 2000 October and were com- 
pleted in 2002 May. The data were taken at the Parkes 64-m 
telescope using the 21 -cm multibeam receiver, which consists 
of an array of 13 beams each with two orthogonal linear po- 
larizations and an average beamwidth of 14'.3 (FWHM). The 
multibeam correlator has a bandwidth of 64 MHz and covers 
the velocity range -1200 < cz < 12700 km s" 1 . The channel 
spacing is 13.2 km s" 1 and the average system temperature is 
20K. 

The northern ZOA region referred to in this paper lies at de- 
clinations > 0° and is divided into four fields, each approxi- 
mately 8° x 10°, as shown in Figure 1 . Each field was scanned 
425 times, with a single 8° scan covering a region 1 .7° wide in 
eight minutes. The central beam for all scans lies along a line 
of constant Galactic latitude, with each of the 425 scans being 
separated by 1 .'4 in latitude. At the midpoint of each scan, the 
feed angle is 15° with respect to the scan direction so that all 

6 The Parkes telescope is part of the Australia Telescope which is funded 
by the Commonwealth of Australia for operation as a National Facility man- 
aged by CSIRO. 



13 beams lie at different latitudes in a minimally redundant 
fashion. The total integration time is 2100 s beam" 1 . The me- 
dian rms noise of the four fields after Hanning smoothing is 
6.0 mJy beam" 1 . 

The data were bandpass-corrected, Doppler-corrected, and 
calibrated using AIPS++ LIVEDATA (Barnes et al. 2001). 
Gridding was performed by GRIDZILLA, using a top-hat me- 
dian gridding algorithm with a top-hat radius of 6' (Barnes et 
al. 2001). The resulting pixel and beam sizes are 4' x 4' and 
15'.5, respectively. Spectral ringing associated with Galac- 
tic HI emission and the baseline ripple associated with strong 
continuum emission were respectively removed or suppressed 
using Hanning smoothing and the 'scaled template method' 
(see Barnes et al. 2001). Hanning smoothing increases the 
velocity resolution to 27 km s" 1 . 

Several sources of extraneous signal remained after stan- 
dard processing. Residual baseline ripples caused by con- 
tinuum radio sources are the most prominent. Hydrogen re- 
combination lines at 1424.734 MHz (H166a), 1399.368 MHz 
(HI 67a), and 1374.601 MHz (HI 68 a) appear at velocities of 
cz=-91 1 km s" 1 , 4507 km s" 1 , and 9990 km s" 1 , respectively. 
Interference from the L3 beacon of the GPS near 1380 MHz 
(8778 km s" 1 ) is also present. These sources of interference 
can generally be identified by their characteristic frequencies 
and widths. In the case of recombination lines, they are ac- 
companied by continuum emission from their underlying HII 
regions. 

3. GALAXY SEARCH 

The four Hanning-smoothed data cubes were visualized us- 
ing KVIEW in the software package KARMA (Gooch 1995) 
and were searched by eye; the strong residual baseline rip- 
ples prevented the use of automatic galaxy-finding software. 
The galaxy search was conducted over the full velocity range, 
-1200 < cz < 12700 km s" 1 , although confusion from Galac- 
tic emission prevented a complete survey of the region with 
| v | < 250 km s" 1 . All possible planes, RA-Dec, RA-velocity, 
and Dec-velocity were examined by three independent parties 
and a final galaxy list was decided upon by a fourth author. 

3.1. Determination of HI Properties 

The coordinates and spectral characteristics of each galaxy 
were measured using the program MBSPECT in the software 
package MIRIAD (Sault, Teuben, & Wright 1995). The galaxy 
positions were determined by spatially fitting a Gaussian to 
the velocity-integrated emission of each galaxy. We subse- 
quently used this position in the spectral analysis, where the 
weighted sum of the emission in each velocity plane was used 
to form the spectral profile. The total flux was measured by 
integrating the flux over the galaxy profile. The systemic ve- 
locity was taken to be the midpoint of the profile at 50% of the 
peak flux. The widths at 20% and 50% of the peak flux density 
were measured using a width-maximizing technique. Due to 
the large Hanning-smoothed velocity resolution of 27 km s" 1 , 
the velocity widths we measure are broader than the intrin- 
sic widths of the galaxies. The velocity widths at the 20% 
and 50% peak flux levels were therefore corrected by 21 and 

14 km s" 1 , respectively (see Henning et al. 2000). If the HI 
emission showed no evidence of extension or confusion, we 
analyzed a region 28' around the galaxy. For seven of the 
galaxies, a narrower box of width 12' or 20' was used to iso- 
late the galaxy emission from that arising from a neighboring 
galaxy, which may have led to an underestimation of the total 
flux. While no galaxies in our survey are clearly extended, we 
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FIG. 1 . — The survey region consists of four fields, each approximately 8° X 



did identify one galaxy that is possibly extended or confused, 
J1912+02. The spectral fitting for this galaxy was conducted 
by summing the emission over a 36' x 36' box. The position 
was fit over a smaller 12' x 12' box. 

Where possible, a first-order baseline was fit to each pro- 
file. For approximately half of the galaxies, however, base- 
line ripples made a first-order fit unrealistic; for these galax- 
ies, we used the lowest order baseline that gave a reasonable 
fit. The distance to each galaxy, D = vlg/H , was calcu- 
lated using its velocity in the Local Group standard of rest, 
vlg = Vhei + 300 sin/ cos/?. The HI mass was determined using 
M ra [M Q ] = 2.36 x 10 5 D 2 F ffl , where D is the distance in Mpc 
and /"hi is the integrated HI flux in Jy km s" 1 . 

To estimate the errors on the HI parameters, we conducted 
a set of simulations. Each galaxy profile was first smoothed 
using a Savitzky-Golay smoothing filter (Press et al. 1992, 
Sec. 14.8). This least-squares polynomial filter reduces noise 
while better retaining the higher-order characteristics of the 
spectrum than traditional smoothing techniques. Fifty sim- 
ulated spectra were then created for each galaxy by adding 
random Poisson noise to the smoothed galaxy spectrum; the 
rms of this random noise was set equal to that of the original 
galaxy spectrum. The errors on the HI parameters were then 
taken to be the median absolute offsets between the param- 
eters of the 50 simulated spectra, as measured by MBSPECT, 
and those of the smoothed galaxy profile. The median errors 
on the peak and integrated flux for the entire sample of galax- 
ies are approximately 9% and 4%, respectively. The median 
error on the systemic velocity is 3 km s -1 and the median er- 
rors on the 50% and 20% velocity widths are 8 km s -1 and 
10 km s" 1 , respectively. We note that the error we quote for 
the integrated flux does not take into account the uncertainty 
due to baseline subtraction or the uncertainty in calibration; 
the total errors on the integrated flux are likely to be ~ 10- 
15%. We further discuss this method and the dependence of 
the estimated errors on S/N in Appendix A. 

3.2. Results 

11 definite or probable galaxies were detected, 35 in the 
/ = 36° to 52° region and 42 in the / = 196° to 212° region. 
Because adjacent data cubes extend slightly beyond the 8° x 
10° boundaries and overlap in Galactic longitude, 5 galaxies 
were initially detected twice. Of these, the detection located 
within the well defined 8° x 10° region of a cube was chosen 
for the final list. The HI properties of the 77 galaxies are 
presented in Table 1 . The HI spectra of a subset of the sample 
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are shown in Figure 2; the complete set of spectra is available 
in the electronic edition of the Journal. The set of dotted lines 
nearest the galaxy spectrum defines the region in which the 
galaxy profile was measured. A second set of more widely 
spaced dotted lines indicate the region over which the baseline 
was fit. If no second set is present, the baseline fit was applied 
to the entire region displayed. All velocities are given in the 
optical convention, v = cz. 

For several galaxies, a highly variable baseline prevents a 
convincing spectrum from being attained although the galaxy 
is clearly present in the data cube. As an example, consider 
J0623+14, a galaxy that was detected despite its superposition 
on one of the most extreme ripples that remain in the survey 
data. Coincidentally, J0623+14 has previously been observed 
by Arecibo in an HI follow-up survey of Weinberger galax- 
ies (Pantoja et al. 1994). Despite the extremely poor base- 
line we observe, our HI parameters are in satisfactory agree- 
ment with the Arecibo HI parameters. The Arecibo flux is 
higher than what we detect by 14% and the peak flux density 
is lower by 22%. The most significant difference between our 
HI properties and those derived from the Arecibo observation 
is the 50% velocity width; we measure a value of 286 km s -1 
whereas Pantoja et al. measure 411 km s" 1 . 

The identification of confused and marginally extended 
galaxies in our sample is difficult. We expect most normal 
galaxies to be unresolved in the Parkes beam at the typical dis- 
tances in this survey, although confusion in the Parkes beam is 
not unlikely. The HI profile is the best indicator of a confused 
signal, although a profile can be intrinsically asymmetric, and 
minor contributions to the HI profile often do not disturb the 
profile greatly. A possible confusion in the detection can be 
confirmed by the relative positions and distances of the poten- 
tial optical and IR counterparts, although background galax- 
ies and Hi-poor counterparts make this difficult, as does the 
high extinction in the ZOA. The unambiguous identification 
of confused galaxies can therefore only be made by high- 
resolution follow-up observations. 

The mass distribution of the survey galaxies is shown in 
Figure 3. The HI masses range from 1.0 x 10 7 M© to 
3.6 x 10 10 M Q , with a mean mass of 4.3 x 10 9 M Q . The 
distributions of velocities in the Local Group (solid line) and 
heliocentric (dotted line) standards of rest are shown in Fig- 
ure 4. By using the Local Group standard of rest, we remove 
the effect of the Milky Way's velocity within the Local Group 
from our velocity and distance determinations. As we are in- 
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FIG. 2. — HI spectra of a subset of the survey galaxies. The spectra of all 77 galaxies are available in the electronic edition of the Journal. The inner dotted 
lines indicate the range in velocity over which the galaxy profile was measured. The outer dotted lines indicate the region over which the baseline was fit. 
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FIG. 3. — Mass distribution of galaxies detected in the survey. The mean 
mass is 4.3 X 10 9 Mq. 



terested in large-scale structures outside of the Local Group, 
we feel that this is the more appropriate standard of rest. Both 
the velocities in the Local Group and heliocentric standards of 
rest have been plotted in Figure 4, however, for comparison. 
Distinct peaks can be seen at Local Group standard of rest ve- 
locities near 3250 km s" 1 , 4750 km s" 1 , and 6750 km s" 1 in 
the I = 36° to 52° region which seem to be due to filaments 
crossing the Galactic plane at those velocities (see also Fig. 
10). The first of these peaks seems to be due to the boundary 
of the Local Void. The velocity distribution for the / = 196° 
to 212° region seems smoother overall, except for the curi- 
ous gap around 4500 km s" 1 . As will be discussed in §5, this 
structure may be due to the wall of the Canis Major Void. 

3.3. Nearby Galaxies 

Of the 77 galaxies detected in this survey, three have sys- 
temic velocities of cz < 700 km s" 1 . The newly detected 
galaxy HIZOA J0630+08 (l,b = 203.1° -0.9°) has a veloc- 
ity of only 367 ± 1km s" 1 (see Table 1) and could be a high 
velocity cloud, though its isolation and relative compactness 
suggest it is more likely a nearby dwarf galaxy. The Galac- 
tic extinction towards this galaxy is Ab ~ 3 mag (Schlegel et 
al. 1998). HIZOA J0630+08 may be part of the Orion Group 
(Giovanelli & Haynes 1981) whose known members are the 
Orion Dwarf (l,b = 200.6°, -12.3°), the irregular galaxy 
0554+07 (l,b = 200.0°, -8.4°), the dwarf irregular galaxy 
CGCG 422-003 (l,b = 198.8°, -12.2°) and UGC 03303 (l,b 
= 198.6, -16.9°). Their velocities are cz = 365, 411, 441 
and 521 km s" 1 , respectively (Giovanelli 1979, Giovanelli & 
Haynes 1981, Michel & Huchra 1988, Thuan & Seitzer 1979). 
Karachentsev & Musella (1996) give distances of 6.4 ±2.2 
Mpc (Orion Dwarf) and 5.5 ± 1.9 Mpc (0554+07) assuming 
Galactic extinctions of 2.7 and 2.9 mag, respectively. At a 
distance of 6 Mpc the HI mass of HIZOA J0630+08 would be 
4 x 10 7 M Q . Further investigations of the Orion Group and 
surroundings are under way. 

In addition to HIZOA J0630+08, we detect two other 
nearby galaxies, HIZOA J1940+11 and HIZOA J1914+10, 
with systemic velocities of 580 km s -1 and 654 km s" 1 , re- 
spectively. For HIZOA J 19 14+ 10 we derive an HI mass of 
M H i = 6.5 x 10 8 M©. This galaxy was previously detected 
in the Dwingeloo Obscured Galaxies Survey and is known as 
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FIG. 4. — Velocity distribution for galaxies in the / = 36° to 52° and 
I = 196° to 212° regions. The velocity distributions in the Local Group and 
heliocentric standards of rest are given by a solid line and a dotted line, re- 
spectively. 



Dw044.8-0.5 (Henning et al. 1998, Rivers 2000). The Galac- 
tic extinction in this direction is extremely high, preventing 
detection of optical or infrared counterparts. 

For HIZOA J 1940+1 1, a newly cataloged galaxy, we derive 
an HI mass of only M H i = 4.1 x 10 7 M . Confirmation of 
this HI detection would be desirable as it is very weak and 
found at the edge of our survey field. The Galactic extinction 
in this direction is Ab ~ 2 mag (Schlegel et al. 1998). Both 
HIZOA J1940+11 and HIZOA J1914+10 appear isolated, in 
contrast to HIZOA J0630+08, and lie in the outskirts of the 
Local Void. 

3.4. Survey Completeness 

Our ability to detect a galaxy in the survey depends on both 
its integrated flux and velocity width. A galaxy with a high 
total flux and a small to moderate velocity width can be iden- 
tified most easily; decreasing a galaxy's flux or increasing its 
velocity width would generally make detection more difficult. 
This, and previous experience with HI surveys (e.g. Kilborn 
et al. 2002), leads us to expect the survey to be approximately 
mean flux limited, where mean flux is defined here to be the 
integrated flux divided by the 20% velocity width (or the 50% 
velocity width in cases where the 20% width is not available). 
This expectation is confirmed by plotting the number of sur- 
vey galaxies with a given flux, N(S), where S is defined to 
be the mean, integrated, or peak flux density. We find that 
N(S) declines most steeply at low fluxes, and therefore best 
defines the flux cutoff, when the flux is taken to be the mean 
flux density. 

A plot of N(S) against S is shown in Fig. 5, with S repre- 
senting the mean flux density as above. It is consistent with a 
Euclidean power law, N(S) oc S~ef n , and a mean flux density 
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FIG. 5. — Number of survey galaxies with a given mean flux density. A 
Euclidean power law, N(S) oc S^^ n , has been fit to the data. The approximate 
mean flux completeness limit, 0.022 Jy, was determined by averaging the fit 
of this power law to regions of low, medium, and high rms noise. 



completeness limit of 0.022 Jy, assuming approximate spa- 
tial homogeneity. The best-fit linear relationship between the 
mean flux density and peak flux density of galaxies in the sur- 
vey has a slope of 0.75. The above completeness limit there- 
fore corresponds roughly to a peak flux density completeness 
limit of 0.029 Jy, which is ~ 5a above the rms noise of the 
survey, 6.0 mJy. 

In Figure 6, we plot as a solid histogram the area of the sur- 
vey, A(< S), at or below a given rms noise, S. The number of 
galaxies detected at or below a given rms noise is also given 
as a dotted histogram. We conservatively define the effective 
area of the survey to be the area with an rms noise at or be- 
low 0.0073 Jy, 239 deg 2 . Inside this area, the limiting mean 
flux density of the survey, 0.022 Jy, corresponds to a 3<r, or 
higher, detection (corresponding to 5<r, or higher, in peak flux 
density). This area is less than the nominal survey area of 
4 x 8° x 10° = 320 deg 2 partly because of edge effects, and 
partly because of the increased continuum (particularly resid- 
ual continuum ripple) along the Galactic Plane. 

To investigate the effect of the variable rms, we divide the 
survey into three regions of equal area, but of different rms 
noise. We again plot N(S) against S for each of the three re- 
gions and fit the differential number counts with a Euclidean 
power law. We estimate completeness limits of 0.019, 0.020 
and 0.026 Jy for the regions of low, medium, and high rms 
noise, respectively. As expected, the completeness limit in- 
creases with rms noise. However, the variation is sufficiently 
limited that, for the present purposes, a mean completeness 
limit of 0.022 Jy adequately characterizes the survey. 

Finally, it is worth noting that, in addition to good com- 
pleteness above a mean flux density of 0.022 Jy, this survey 
is expected to be highly reliable. Galaxy identifications were 
made by three authors independently, and the final galaxy list 
was made by a fourth author in a conservative manner, guided 
by experience gained from follow-up observations of surveys 
with similar noise characteristics (Henning et al. 2000, Zwaan 
et al. 2004). We expect that very few, if any, of the new galax- 
ies presented here will be false detections, and therefore the 
reliability of the catalog to be close to 100%. Ultimately, how- 
ever, follow-up HI observations of some weaker objects will 
be required to ensure that no false detections remain. 

4. MULTIWAVELENGTH COUNTERPARTS 

The region of sky covered by the northern ZOA survey 
overlaps slightly with a number of other surveys. These in- 
clude four blind HI surveys, the Arecibo Dual Beam Sur- 
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FIG. 6. — Histogram of the survey area (solid line, left axis) and number of 
galaxies (dotted line, right axis) detected at or below a given rms noise. The 
survey completeness is set to be the area of the survey with an rms noise at or 
below 0.0073 Jy. 



vey (ADBS; Rosenberg & Schneider 2000), the pilot NRAO 
300-ft telescope survey (Henning 1992), the Dwingeloo Ob- 
scured Galaxy Survey (DOGS; Henning et al. 1998, Rivers 
2000), and the HI Parkes Zone of Avoidance Shallow Survey 
(HIZSS; Henning et al. 2000), although the overlap with the 
HIZSS is extremely small. This region also coincides with a 
number of HI follow-up surveys and several optical surveys, 
from which the ZOAG (see Lercher, Kerber, & Weinberger 
1996), CGMW (see Roman, Iwata, & Saito 2000), WEIN 
(Weinberger 1980), and additional catalogs have been com- 
piled. Significant infrared coverage of these regions has been 
provided by IRAS and 2MASS. 

4.1. Counterpart Search 

We searched for HI, optical, and infrared counterparts to 
the northern ZOA galaxies using two methods. We first used 
the NASA/IPAC Extragalactic Database (NED 7 ), the 2MASS 
Extended Source Catalog (XSC 8 ), and the IRAS Point Source 
Catalog (PSC) to conduct an automated search for the near- 
est potential counterparts to each HI galaxy. The fraction of 
these counterparts likely to be associated with the HI galax- 
ies was estimated by a simulation, described in Appendix B. 
Our second approach was to determine the most likely coun- 
terparts by searching for and visually inspecting all potential 
optical and infrared counterparts, including (1) counterparts 
located further from the HI position than those chosen by our 
previous search and (2) previously unidentified galaxies in the 
Digitized Sky Survey (DSS 9 ) POSS-II near-infrared, red, and 

7 This research has made use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion Laboratory, California Insti- 
tute of Technology, under contract with the National Aeronautics and Space 
Administration. 

8 This publication makes use of data products from the Two Micron All 
Sky Survey, which is a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Center/California Institute of Technol- 
ogy, funded by the National Aeronautics and Space Administration and the 
National Science Foundation. 

9 The Digitized Sky Surveys were produced at the Space Telescope Sci- 
ence Institute under U.S. Government grant NAG W-2166. The images of 
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FIG. 7. — Histograms of the offset between the HI velocities and fluxes measured by our survey, and those from the literature. The mean offset in velocity is 
6 km s -1 and that in flux is -16%. 



blue images. Any galaxy visible in the DSS bands was clas- 
sified by us as an optical detection. The visual search is, of 
course, subject to error; we have tried to minimize this error 
by having multiple authors participate in the choice of coun- 
terparts. Previous investigations of the offset between Parkes 
HI sources and their counterparts (see Juraszek et al. 2000, 
Koribalski et al. 2004, Ryan-Weber et al. 2002) indicate that 
while most counterparts are found within 3', offsets as high 
as 5' are not highly unusual. A conservative search radius of 
10' was used for the initial searches to ensure that no potential 
counterparts were missed. 

4.2. Previously Detected HI Galaxies 

Of the 77 galaxies detected here, 20 have been detected pre- 
viously in HI. We have detected all galaxies from the ADBS, 
NRAO 300-ft, HIZSS, and DOGS blind HI surveys that lie 
in our survey region, although J0635+11 was resolved into 
two galaxies separated by 4' by the ADBS. We have also de- 
tected several galaxies from follow-up HI surveys. Figure 7 
illustrates the difference between our measured systemic HI 
velocities and the HI velocities given by NED as well as the 
fractional differences in flux between our survey and measure- 
ments in the literature (Rosenberg & Schneider 2000, Pantoja 
et al. 1997, Pantoja et al. 1994, Theureau et al. 1998, Hen- 
ning et al. 2000, Rivers 2000). The distribution of velocity 
offsets has a mean of 6 km s" 1 and a dispersion of 19 km s -1 . 
The flux we measure is offset by an average of -16% from the 
values in the literature, with a dispersion of 32%. 

The galaxies with the largest positive and negative flux 
offsets are HIZOA J0654+08 and HIZOA J0622+11B, re- 
spectively. Both were previously observed as part of the 
Arecibo Dual Beam Survey. The small ADBS flux for HI- 
ZOA J0622+1 IB can potentially be explained by the fact that 
ADBS flux measurements were underestimated if the source 
was offset from the beam center (see Rosenberg and Schnei- 
der 2000). We are confident in our fits to both galaxies, and 
are unable to explain the large positive flux offset of HIZOA 
J0654+08. 

Two of the remaining galaxies for which our flux is sig- 

these surveys are based on photographic data obtained using the Oschin 
Schmidt Telescope on Palomar Mountain and the UK Schmidt Telescope. 
The plates were processed into the present compressed digital form with the 
permission of these institutions. 



nificantly higher than that in the literature were identified by 
the Dwingeloo Obscured Galaxies Survey. Both have strong 
clean profiles. The discrepancies in flux most likely arise 
from the conservative way in which the Dwingeloo galaxy 
fluxes were measured (Rivers (2000) defined the total inte- 
grated flux to be the flux inside the N H i = 1.25 x 10 20 cm" 2 
contour). Also, as the Dwingeloo fluxes were measured from 
VLA maps, any extended emission would have been missed 
by the interferometer. 

4.3. Results of Counterpart Search 

We present in Table 2 the HI counterparts from other sur- 
veys of the ZOA galaxies, along with the most likely opti- 
cal, 2MASS, and IRAS counterparts. For the latter group, 
we include only those sources that were visually verified to 
be likely counterparts. Potential counterparts were discarded 
from this list if they did not appear to be a galaxy, were off- 
set from the HI position by a statistically unlikely amount 
(> 4'), or were of an unlikely galaxy type. In seven cases, the 
2MASS counterpart that was chosen by the visual inspection 
was not the counterpart closest to the HI position. All likely 
optical and IRAS counterparts (excluding the DSS galaxies 
which were not found via the automated NED search) were 
the counterparts nearest the HI position. In Table 3, we give 
the positions of the newly detected DSS counterparts. 

The spatial distribution of sources with multiwavelength 
counterparts is shown in Figure 8, where solid circles repre- 
sent galaxies with both optical and infrared counterparts, solid 
triangles represent galaxies with optical counterparts only, 
and solid stars represent galaxies with infrared counterparts 
only. HI galaxies with no counterparts are given as crossed 
open circles. We have overlaid contours of E(B-V) on Figure 
8 using the DIRBE/IRAS data. Contour values of l m .O and 
3™0 were used for the I = 36° to 52° region, whereas values 
of m .3 and m .8 were used for the contours of the / = 196° to 
212° region. 

In summary, 19 of our 77 galaxies have likely optical 
counterparts, 27 have likely 2MASS counterparts, and 11 
have likely IRAS counterparts. In addition, 16 galaxies have 
likely optical DSS counterparts that are previously uncata- 
logued. The percentage of counterparts found via the auto- 
mated search that were still considered likely after a visual 
verification is roughly consistent with the expected percent- 
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FIG. 8. — Spatial distribution of galaxies. HI galaxies with no counterparts are represented by crossed open circles. Solid circles represent galaxies with both 
optical and infrared counterparts, solid triangles represent galaxies with optical counterparts only, and solid stars represent galaxies with infrared counterparts 
only. The counterparts used in this plot are only those that remain after the visual verification. The contours overlaid on the plot are contours of E(B-V) (Schlegel 
et al. 1998), with values of l m .O and 3 m .O for the / = 36° to 52° region, and m .3 and m 8 for the I = 196° to 212° region. 



ages calculated by means of the simulations discussed in Ap- 
pendix B. 

5. LARGE-SCALE STRUCTURE 

The spatial distribution of galaxies in our survey can be 
found in Figures 9 through 12. In Figures 9 and 11, the HI 
galaxies detected in our survey are given by solid points. We 
have plotted as open points all galaxies in the Lyon-Meudon 
Extragalactic Database (LEDA) 10 with known velocities. In 
each velocity-slice, the galaxies with the lowest velocities are 
represented by triangles; medium and high velocity galaxies 
are plotted as squares and circles, respectively. The dotted 
box outlines the region sampled by our survey. In the cone 
diagrams of Figures 10 and 12, the HI galaxies are plotted as 
solid circles and LEDA galaxies are plotted as open circles. 
The survey region is again defined by dotted lines. The veloc- 
ities in Figures 9, 10, 11, and 12 have been converted to the 
Local Group standard of rest. Void parameters are taken from 
Fairall (1998) and are given in Table 4. 

The I = 36° to 52° region lies near several voids, including 
the Local, Microscopium, and Cygnus Voids. At first glance 
it appears that the galaxies detected in our survey may define 
the edge of the Local and/or Microscopium Voids, but we cau- 
tion that preliminary results from the full ZOA survey indicate 
that the number of galaxies rises slightly as one moves from 
our survey region towards the Local and Microscopium Voids 
before dropping off at Galactic longitudes of 10° < / < 20°. 

When visualized including even larger areas (not shown 
here), the nearby detected galaxies that produce the distinct 
peak around 3000 km s -1 appear to be the smooth continua- 
tion of the sine-wave like feature that can be traced across the 
whole southern sky (Kraan-Korteweg, Koribalski, & Juraszek 
1999). This wave feature crosses the Milky Way in Puppis and 
continues towards the Antlia, Hydra, and Centaurus structures 
before folding back across the Galactic Plane in a structure 
called the Centaurus Wall by Fairall (1998). This wave encir- 
cles the Local Void while crossing the ZOA for a third time 
towards the region inspected here. 

The second peak in the velocity distribution around 

10 http://leda.univ-lyonl.fr - the Lyon-Meudon Extragalactic Database is 
supplied by the LEDA team at the Centre de Recherche Astronomique de 
Lyon. 



TABLE 3 
New DSS Counterparts 


HIZOA 


"2000 


^2000 


J0622+11A 


06 21 59.4 


+11 18 38 


J0634+09 


06 34 04.9 


+09 12 55 


J0635+11 


06 35 54.6 


+11 08 10 


J0641+01 a 


06 41 04.4 


+01 50 24 


J0659+06 


06 59 38.1 


+06 27 19 


J0702+03 


07 02 50.4 


+03 11 14 


J1853+09 3 


18 53 47.5 


+09 51 16 


J1917+07 


19 17 24.8 


+07 49 07 


J1927+09 a 


19 27 49.5 


+09 27 32 


J 1929+08 


19 29 19.9 


+08 02 42 


J1930+12 a 


19 30 32.2 


+12 11 41 



"We see a galaxy at this position, but cannot determine whether it 
is the correct counterpart to the HI galaxy. 



TABLE 4 

Void Parameters from Fairall (1998) 



Name 


/ 

(deg) 


b 

(deg) 


cz 

(kms" 1 ) 


radius 
(kms- 1 ) 


Aquarius Void 


60 


-41 


4500 


1500 


Canis Major Void 


229 


-13 


5000 


2500 


Capricornus Void 


2 


-17 


8500 


~2500 


Cor Bor Void 


58 


37 


5200 


2000 


Cygnus Void 


67 


-9 


3500 


1250 


Delphinus Void 


59 


-6 


2500 


1500 


Gemini Void 


172 


9 


3000 


~ 1250 


Local Void 


18 


6 


1500 


1500 


Microscopium Void 


10 


1 


4500 


1750 


Orion Void 


206 


-2 


1500 


750 


Taurus Void 


167 


-29 


4000 


~ 2000 



4500 km s also seems to be due to a clear filament present 
in both the middle panel of Figure 9 as well as the redshift 
cones of Figure 10. Here, the detected galaxies seem to di- 
vide the Microscopium and Cygnus/Delphinus Voids, the lat- 
ter two of which are likely to form one void rather than two 
distinct voids. 
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FIG. 9. — Spatial distribution of galaxies in the / = 36° to 52° region. The 
velocities are in the Local Group standard of rest. The dotted line outlines the 
survey region. HI galaxies detected by this survey are represented by solid 
points; LEDA galaxies are given by open points. In each velocity range, 
triangles, squares, and circles represent the galaxies in the low, middle, and 
high velocity bins, respectively. The dotted line outlines the survey region. 



FIG. 10. — Cone diagrams of galaxies in the / = 36° to 52° region. The 
velocities are in the Local Group standard of rest. HI galaxies detected by 
this survey are represented by solid circles; LEDA galaxies are given by open 
circles. The survey region is defined by dotted lines. 



There are several large-scale structures in the vicinity of 
the I = 196° to 212° region, which lies between the Perseus- 
Pisces supercluster and the Puppis filament. Two of the near- 
est newly-detected galaxies, HIZOA J0700+01 and HIZOA 
J0705+02, are found within 4.8 and 4.9 Mpc of the center of 
the Orion Void. This void (l,b,cz = 206°, -2.0°, 1500 km s" 1 ) 
was previously thought to have a radius of 750 km s -1 , or 
10 Mpc (Fairall 1998). The LEDA Database, however, con- 
tains 42 galaxies that lie between 5 Mpc and 10 Mpc of the 
void center, and two that lie at smaller radii (4.6 and 4.9 Mpc). 
It therefore appears that the true radius of this void is approx- 
imately 5.0 Mpc and that the two galaxies detected here lie 
along its edge. 

The galaxies between 2500 and 5500 km s -1 seem to pro- 
vide the border of two voids, the Gemini Void, now quite 
prominent in the upper panel of Figure 11, and the Ca- 
ms Major Void, which now seems slightly smaller and cen- 
tered at higher latitudes (l,b,cz = 220°, 0°, 5000 km s" 1 ) as 
compared to its location as given by Fairall 1998 (l,b,cz = 
229°, -13°, 5000 km s" 1 ). The distinct gap between 4000 and 
5000 km s" 1 in the velocity distribution (Fig. 4) is due to the 
fact that the continuous part of the border of the Canis Major 
Void between 4000 and 5000 km s" 1 falls just outside of our 
survey regions and crosses the ZOA at quite constant veloci- 



ties. 

6. MASS FUNCTION 

Over the past several years, a number of attempts have 
been made to better define the HI mass function. Most re- 
cently, the HIPASS Bright Galaxy Catalog (BGC, Zwaan et 
al. 2003) and the ADBS (Rosenberg & Schneider 2002) have 
been utilized. Prior to these surveys, which cataloged 1000 
and 265 galaxies, respectively, the HIPASS SCC (Kilborn 
2000), the HIPASS HIZSS (Henning et al. 2000), and the 
AHISS (Zwaan et al. 1997) provided the best estimates to the 
mass function. The sensitivity and coverage of our survey do 
not allow us to estimate the mass function any more precisely 
than these recent measurements, especially at the low-mass 
end where the slope of the mass function is most uncertain. 
Determining the mass function is still of interest, however, 
as it allows an investigation of the comparative mass density 
and distribution of galaxies in our survey region, as well as a 
rough comparison with previous estimations. 

We calculate the HI mass function in two ways, the stan- 
dard 1/Vmax method (Schmidt 1968) and a refined 1/V ma x 
method (Saunders et al. 1990). The latter incorporates the 
density profile of the galaxies into the calculation to reduce 
the effect of clustering and large-scale structure. The stan- 
dard 1 /y max method non-parametrically determines the mass 
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TABLE 5 
SCHECHTER PARAMETERS 



Method 


log 6* 


a 


log M*, 


log HI density 3 


log HI density b 


1/V ma x 

Modified 1 /Vmax 


-2.33 ±0.33 
-2.35 ±0.32 


1.13±0.36 
1.16 ±0.27 


9.85±0.19 
10.05 ±0.20 


7.56 ±0.42 
7.75 ±0.37 


7 <;<:+O.07 



"HI mass density calculated from Schechter parameters (log Mq Mpc 3 ) 
b HI mass density calculated from E^ 1 (log M Mpc -3 ) 



°..5. Y. <. 3600_ (km s-') 




240 220 200 180 160 
Galactic Longitude 



FIG. 1 1 . — Spatial distribution of galaxies in the / = 196° to 212° region. 
The velocities are in the Local Group standard of rest. The dotted line outlines 
the survey region. HI galaxies detected by this survey are represented by solid 
points; LEDA galaxies are given by open points. In each velocity range, 
triangles, squares, and circles represent the galaxies in the low, middle, and 
high velocity bins, respectively. The dotted line outlines the survey region. 




FIG. 12. — Cone diagrams of galaxies in the I = 196° to 212° region. The 
velocities are in the Local Group standard of rest. HI galaxies detected by 
this survey are represented by solid circles; LEDA galaxies are given by open 
circles. The survey region is defined by dotted lines. 



function by weighting each galaxy's mass by the inverse of 
the maximum volume in which it could be detected by the 
survey. This method, however, is potentially affected by 
large-scale structure, as it assumes that the galaxy popula- 
tion has a homogeneous distribution. The Saunders et al. 
method is able to remove the majority of the effects of large- 
scale structure by first using a maximum likelihood technique 
to determine the radial density distribution of the sample, 
p(r). This density profile is used to define an effective vol- 
ume, V e ff = Jo™"' 1 r 2 p(r)dr, which replaces V max in the 1 /V m . dx 
method. In both cases, the resulting mass function can be fit 
by a Schechter (1976) function: 



0(M m )dM m = 9*(-^Texp(--^)dM H1 

The HI mass function can be found in Figure 13, where the 
data points from the standard 1 /V max method are given by cir- 
cles and the fit by a solid line. The modified 1 /V max method is 
represented by triangular points and a dashed line. For com- 
parison, the HI mass function determined from the BGC has 
been included on the plot, as a dotted line. Only those galax- 
ies with a mean flux above our completeness limit of 0.022 Jy 
were included in the mass function calculation. A bin size of 
0.5 dex was used for both techniques. For the standard 1 /V max 
method, the points are plotted at the mean mass of the galaxies 
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FIG. 13. — (a) HI mass function of the survey galaxies. The circles represent the mass function calculated from the standard 1 /V ma x method. The solid line is 
a least-squares fit of a Schechter function for this method. The modified 1 /V max method is represented by triangular points. The Schechter function fit for this 
method is given by a dashed line, (b) HI mass function of the two survey regions. The mass function of the / = 36° to 52° region is represented by circular points 
and a solid line. The mass function of the / = 196° to 212° region is represented by triangular points and a dashed line. The bin size is 0.5 dex; bins are plotted 
if one or more galaxies are present. Points are plotted at the mean mass within the bin for the standard method, and at the mean bin position for the modified 
method. For both plots, the dotted line represents the mass function from the HIPASS BGC (Zwaan et al. 2003). 



in each bin. Missing points indicate that we detect no galax- 
ies in that bin. For the modified 1 /V max , the points represent 
the middle of the mass bin. The galaxy J0635+1 1, which was 
resolved into two galaxies by the ADBS, does not affect our 
mass function, as its mean mass falls below our cutoff. For a 
discussion on the effects of confusion on HI mass functions, 
see Zwaan et al. (2003). 

The Schechter parameters from the two methods as well 
as the derived HI mass densities can be found in Table 5. 
The two methods give similar mass functions, suggesting that 
large-scale structure in these regions does not strongly affect 
the calculation of the mass function. The largest deviation 
between the two methods occurs at the high-mass end, where 
the standard method is in better agreement with the BGC mass 
function than the modified method. In addition, the standard 
1 /V max method gives a < y— > of 0.54 ±0.04, indicating that 
spatial inhomogeneity will not affect the results greatly. Both 
methods of estimating the mass function give a function that 
falls below that of the BGC at the low mass end, with the 
largest discrepancy coming from an underdensity of galaxies 
of mass M ~ 10 8 - 10 9 M Q . This offset may be due to an 
underdensity of mass in the regions examined. Zwaan et al. 
(2003) have shown that the four quadrants of the BGC give 
mass functions with similar values of M* but with varying 
slopes and normalizations with a degree of variation similar 
to what we observe between our mass function and that de- 
rived from the BGC. 

We also include in Figure 13 a plot of the mass functions 
for the two separate survey regions. The I = 36° to 52° re- 
gion is given by circular points and a solid line. The I = 196° 
to 212° region is given by triangles and a dashed line. Due 
to the low number of galaxies, these mass functions are not 
highly reliable, but can still be used for general comparison 
purposes. The mass function for the / = 196° to 212° region 
fits the BGC mass function remarkably well in the high-mass 
regime, whereas the / = 36° to 52° region mass function pro- 
vides the best fit to the BGC mass function in the low mass 
regime. 

7. SUMMARY 



We have detected 77 HI galaxies in the northern extension 
of the Parkes ZOA survey. The survey has a median rms noise 
of 6.0 mJy beam" 1 , is approximately complete to a mean flux 
density of 22 mJy, and has an effective area of 239 deg 2 . Of 
the 77 galaxies detected here, 19, 27, and 11 have likely op- 
tical, 2MASS, and IRAS counterparts already cataloged, re- 
spectively. In addition, a further 16 have likely optical coun- 
terparts visible in the Digitized Sky Survey. Twenty have 
been previously detected in HI. The spatial and velocity dis- 
tributions of the galaxies in the I = 36° to 52° region reveal 
several filaments that cross the Galactic plane, the nearest of 
which seems to be the continuation of a sine-wave like fea- 
ture that can be traced across the whole southern sky. The 
medium-velocity galaxies in the / = 196° to 212° region seem 
to provide the border of the Gemini and Canis Major Voids, 
whereas the detection of two new low-velocity galaxies within 
the Orion Void help to redefine its radius. Of particular inter- 
est is the galaxy HIZOA J0630+08 (l,b = 203°, -0.9°), which 
has a velocity of 367 ± 1 km s -1 in the heliocentric standard 
of rest. We suggest that it belongs to the nearby Orion Group 
which includes a small number of dwarf galaxies. An HI mass 
function was derived using both the standard 1 /V max method 
and a revised 1 /V max method that is less sensitive to large- 
scale structure. Both methods produce similar mass functions 
whose normalizations fall slightly below that of the HI Parkes 
Bright Galaxy Catalog. The observations for the southern 
Parkes ZOA survey are complete and analysis is currently un- 
derway. Upon completion, this survey will be combined with 
the northern extension studied here, which should consider- 
ably increase our understanding of the structure and dynamics 
of the obscured Local Universe. 
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and CONACYT research grant 40094-F (RCKK and JMII). 
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APPENDIX 

A. ERRORS IN HI PROFILE PARAMETERS 

To estimate the errors on the HI parameters, each galaxy profile was smoothed using a Savitzky-Golay smoothing filter with 
a smoothing width of 8 velocity channels (Press et al. 1992, Sec. 14.8). Fifty simulated spectra were then created for each 
galaxy by adding random Poisson noise to the smoothed galaxy spectrum; the rms of this random noise was set equal to that of 
the original galaxy spectrum. The errors on the HI parameters were then taken to be the median absolute offsets between the 
parameters of the 50 simulated spectra and those of the smoothed galaxy spectrum. 

To better understand the errors on the HI parameters and the error simulation technique, we investigate below three additional 
aspects of the error estimation: (1) the effect of S/N on the measured errors, (2) the accuracy of the smoothing filter in preserving 
the galaxy parameters and (3) the biases introduced to the errors by the simulation process. We explore these issues using 
the spectrum of J 1853+09, which has a S/N of 24. This spectrum, and the smoothed profile created using the Savitzky-Golay 
smoothing filter, are shown in Figure Al. A spectrum with a high S/N was chosen to ensure that the smoothed spectrum, used 
below as the template, was as representative as possible of a real galaxy profile. 

To explore the effect of S/N on the errors, we began by running the error simulation, as described in §3.1, on the template. 
The simulation was run 21 times; on each occasion, the rms of the random Poisson noise was adjusted to give a desired S/N. The 
resulting set of 21 simulations have S/N ratios ranging incrementally from 3 to 24. The errors on the HI parameters, as a function 
of S/N, are given in Figure A2(a), where the points represent the median errors and the error bars represent the median absolute 
value of the errors (in both the +y and -y directions). The errors begin to rise at a S/N of 10, and do so quite substantially for S/N 
< 5. This result is in agreement with that of Roth, Mould, & Staveley-Smith (1994), who performed a similar set of simulations 
and found that the errors on HI line widths are particularly unreliable for S/N < 5. Of the 77 galaxies in our survey, none have a 
peak S/N < 5, 21 have 5 < S/N < 10, and 56 have S/N > 10. 

Each of the 21 error simulations described above created 50 simulated galaxy profiles with a given S/N. To investigate the 
ability of the smoothing filter to preserve the galaxy profile, we next smoothed each set of 50 simulated profiles using the 
Savitzky-Golay smoothing filter. The HI parameters of the smoothed spectra were then measured by MBSPECT. Figure A2(b) 
illustrates the difference between the HI parameters of the spectra for which noise was added and then removed and the HI 
parameters of the original smoothed profile. There are several systematic offsets. For S/N > 6, the peak flux density from the 
simulated profiles is lower than the peak flux density of the original smoothed profile, as expected from the smoothing process. 
This offset is quite low, however, tending towards 5% at high S/N. The total flux, systemic velocity, 50% velocity width, and 
20% velocity widths are all slightly larger in the simulated smoothed profiles, but remain relatively constant down to S/N ratios 
of 5, below which the differences in peak flux density, total flux, and 20% velocity width increase. While it is not surprising that 
the smoothing filter increases the total flux and velocity widths, the increase in the systemic velocity is unexpected. This change, 
however, is due to the particular spectrum being used here, and will likely vary with different profiles. As can be seen in Figure 
Al, the overestimation of the velocity at the high-velocity end of the profile is larger than the underestimation at the low-velocity 
end, leading to a slight (<~ 2 km s" 1 ) overestimation of the systemic velocity. 

We lastly investigate how the changes in smoothed parameters discussed above affect the error estimation. To do so, the 
error simulations were run on each of the 50 simulated galaxy profiles with a given S/N ratio. Recall that these 50 simulated 
spectra were created by adding random noise with a given rms to the template. By running the simulations on these 50 profiles, we 
remove the noise that we added (the step explored in the last paragraph) and then re-add it to the smoothed spectra. If no biases are 
introduced by the simulation process, the errors now measured should be equal to the errors estimated by running the simulations 
on the original spectrum of J 1853+09. Figure A2(c) illustrates the differences between the errors derived from the simulated 
profiles and those estimated from the template. The largest bias introduced by the simulation process is the underestimation of 
the error on the 20% velocity width for spectra with low S/N ratios. At S/N > 10, the error on the peak flux is overestimated by 
~ 1.6% and the error on the total flux is underestimated by ~ 0.8%. The errors on the systemic velocity, 50% velocity width, and 
20% velocity width are offset by an average of -0.2 km s" 1 , -0.9 km s" 1 , and 0.2 km s" 1 . At moderate to high S/N, these offsets 
are all well within the median errors measured for the galaxies in this sample (see §3.1), confirming the reliability of this process 
in estimating errors. 



B. COUNTERPART LIKELIHOOD SIMULATIONS 

To investigate the likelihood that the nearest potential optical and infrared counterparts within both 3' and 5' of the HI position 
are real counterparts of the HI galaxies, a simulation was performed. In this simulation, fiducial HI detections were placed across 
our survey region. The positions of these simulated detections were defined by a grid of 120 positions at 20 arcminutes intervals 
surrounding each of the 77 galaxies. Using NED, the 2MASS XSC, and the IRAS PSC, counterparts of these false objects 
were searched for. By measuring the number of objects near these random positions, we determined the number of spurious 
counterparts expected for our real galaxies, and thus the probability, P, that the potential counterparts we detect are in fact related 
to the HI galaxies. 

p _ /V counterparts — N unrelated (B 1 ) 

/Vcounterparts 

_ /V s i mu l a ted galaxies with a potential counterpart ,„,, 
^unrelated — ~ JVgalaxies in survey region \"^) 

^simulated galaxies 

We note that galaxies with low fluxes or large extents often have higher positional errors than bright and/or compact galax- 
ies; the positional uncertainty of a given source is roughly the gridded Parkes beam, 15.5', divided by the signal-to-noise. The 
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FIG. Al . — Spectrum (dotted line) and Savitzky-Golay smoothed profile (solid line) of J1853+09. 
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FIG. A2. — (a) The errors on measured HI parameters for the template as a function of S/N; (b) the difference between the parameters of the smoothed simulated 
profiles and the original parameters of the template; (c) the difference between the absolute value of the HI errors measured from the simulated spectra and those 
measured from the template. The points represent the median errors and the error bars (in both the +y and -y directions) represent the median absolute value of 
the errors. 
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TABLE Bl 

Results of Automated Search for Multiwavelength Counterparts 
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Optical 


OA/f ACQ 


ip AC 
IKAo 


1 = 36° to52°, Ar<5': 








Ncounterparts found via automated search 


2 


6 


24 


% expected to be real 


5% 


77% 


18% 


Ncounterparts that remain after visual inspection" 


1(50%) 


4(67%) 


2(8%) 


/ = 36°to52°,Ar<3': 








Ncounterparts found via automated search with Ar < 3' 


1 


5 


14 


% expected to be real 


26% 


88% 


35% 


Ncounterparts that remain after visual inspection" 


1(100%) 


3(60%) 


2(14%) 


/= 196° to 212°, Ar < 5': 








Ncounterparts found via automated search 


17 


31 


15 


% expected to be real 


96% 


59% 


52% 


Ncounterparts that remain after visual inspection" 


17(100%) 


16(52%) 


9(60%) 


1 = 196° to 212°, Ar < 3': 








Ncounterparts found via automated search with Ar < 3' 


14 


26 


13 


% expected to be real 


96% 


79% 


78% 


Ncounterparts that remain after visual inspection" 


14(100%) 


14(54%) 


9(69%) 



"#(%) gives the number and percentage of counterparts found via the automated search that remain on our list of likely counterparts (and are therefore in Table 3) following a visual 
inspection of all possible counterparts. In seven cases, the most likely 2MASS counterpart was not the nearest 2MASS counterpart to the HI position, and therefore was not the one 
found during our automated search. DSS galaxies were also not found via the automated search, nor was PGC 3097165. 



probabilities given by these simulations are for the entire sample; the above effect will cause slight variations among individual 
galaxies. The simulations were run separately for the two regions of the survey as the density and therefore reliability of counter- 
parts is field-dependent. The difference between simulations with different degrees of clustering around the real HI sample was 
used to evaluate the approximate errors on the estimated counterpart likelihoods, P, <~ ±10%. 

We present in Table Bl the results of the counterpart likelihood simulations. The percentage of counterparts found via the 
automated search that were still considered likely after a visual verification is roughly consistent with the expected percentages 
calculated by means of the simulations. The number of likely counterparts to the galaxies in the I = 196° to 212° region is 
substantially higher than that of the galaxies in the I = 36° to 52° region. The difference between the number of potential 
counterparts in the two survey regions can be attributed to the lower extinction in the latter field. Using the 100 ^m DIRBE/IRAS 
data (Schlegel et al. 1998), we calculate mean B-band extinctions of 12 m .O and 3 m .2 over the / = 36° to 52° and / = 196° to 212° 
regions, respectively. This point is further illustrated by the E(B-V) contours plotted on Figure 8. Although this data are not well 
calibrated in the Galactic plane, they are still useful for general comparison purposes. 
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TAB LEI — Continued 
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The error we estimate for the integrated (lux does not lake into account the uncertainty clue to baseline subtraction or the uncertainty in calibration; the total error on the integrated flux is likely to be ~ 10-15% 
I be 11 ri x v. as calculated over a smaller region to isolate the galaxy emission l'roni that arising from a nearby galaxy. This may have led to an underestimation of the lolal flux. 
The 20% velocity widih could not be accurately measured. 
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from DIRBE//R45 100 /im exlinelion maps (Schlcgcl el al. 1998) 
k Ai' = v' C0lin[er p iirl - v (km s" 1 ) 

C The references for conn terparts are as follows: (ADBS) Rosenberg & Schneider 2000; (CAP) Pantoja et al. 1997; (Dw) Henning el al. 1998. Rivers 2000; (HIZSS) Henning ci al. 2000; (PGC) http://leda.univ-lyonl .fr/ (UGC) Nilson 1973, (WEIN) Weinberger 1980; (ZOAG) Seeberger, Saurer. & Weinberger 1996; (2MASX 
J19061820+1256195)Jarrettetal. 2000 

^ There are two possible counterparts to this source. 2MASS J06254 183+0346338 at a separation of 23' and a DSS galaxy at a = 06 h 25 m 45 s l. 5 = 03 o 47'55" (separation = 1.4'). We do not have enough information to decide between the two. 
e We note thai there is another possible 2MASS counterpart. J06303033+0240026, at a separation of 2.7'. This candidate appears to be of earlier type than that chosen above, but can not be excluded, 
f , 

Our HI detection was resolved into iwo gLikixics squirmed by 4 by the ADBS. We are unable to say to which HI galaxy ilie cuuiUci-parl belongs. 

^ In addition to the counterpart listed here, we see another DSS galaxy al a = 06 h 59 m 38 f i] and S = 06° 25' 19" as well as lwo2MASS candidates: J06595034+0623322(sep=4.3') and J06595034+0623322 (sep=4.4'). From morphology and separation arguments, we believe the candidate we have chosen is the most likely, but we 
can not rule out those listed here. 



